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NOMENCLATURE 


C = Constant 

D = Flexural rigidity Et^/I2(l - V^) 

E = Elastic modulus 

F = Shear force at x = 0 

h = Film thickness 

ho = Value of h where (dp/dx) = 0 

h = h/(6 + 5 ) 

s o 

= (h - ho)/ho 
i = 

= Diffusion length, /Rt/ [ 12( 1-V^ ) ] 

Lj. = X extent of seal lift-off 
Lr = -(Lr/Jlg) 
p = Pressure 
Pg = Sealed gas pressure 

Ps = Et6g/R^ 

Po = Oil pressure; also p at start of film 
P = - Po)/(Ps Pg ■ Po) 
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Po ~ (pg ~ Po^/po 

Rs = Radius of shaft 

Ri = Inner radius of seal 

R = Effective radius, Ri + (t/2) 

t = Thickness of seal wall 

U = Shaft linear velocity 

X = Horizontal coordinate (in direction of U) 

Xm = Location of maximum pressure 

X = -(x/H) 

e 

w = Displacement coordinate (normal to U) 

w = - (w - 6 )/(6 + 6 ) 

o s o 

6m U 6 6y U R^ 

B = 5 

ho^ Pg Et h^^3 

a _ 6m U £ R^ 

'^l e 

Et (6 6 ) 

o s 

6^ = Approximate value of B * 

6 = Interference (positive or negative) i,. 

6^ = Interference due to (pg - Po)» R^(pg ~ Po)/Et 

6 = Interference due to shaft 

s 
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r = ring 
j = index 
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SUMMARY 


An analysis has been developed to predict the performance of pumping Lenin- 
grader reciprocating rod seals when the inlet (gas) side of the seal bore is 
formed by an expansion ring rather than by machining. The prediction of seal 
performance is based on the use of charts which provide necessary design 
parameters without the need for computerized calculations. A numerical exam- 
ple has been included to demonstrate the use of the design charts. Potential 
means for controlling and optimizing both performance and life capability is 
also provided. 

An experimental study was also conducted as part of the overall program in 
order to evaluate pumping Leningrader seals constructed in accordance with the 
analysis. Several seal materials, ranging in elastic modulus from 1.39 x 10^ 
MPa for filled PTFE to 4.58 x 10^ MPa for poly(amide-imide) were tested. For 
the limited number and durational tests conducted, the experiments s.iowed that 
the ring expanded inlet type seal has the ability to provide desired levels of 
sealing. Additionally, it was shown that the lower modulus materials are more 
easily fabricated into seals and provide better sealing capability. 



1.0 INTRODUCTION 


The Pumping Leningrader Seal has been used successfully as a rod seal for . 

Automotive Stirling Engine (ASE) applications. It is used for sealing the 
working fluid (hydrogen) at pressures to 15 MPa against an ambient lubricating 
oil. An analysis was performed to obtain quantitative predictions of the 
mechanism by which the seal operates and to suggest potential means of 
controlling and optimizing seal performance and life. 

The analytical techniques divide the seal into three regions: an inlet zone, a 
contact zone, and an exit zone. The film thickness is largely determined by 
inlet zone behavior, and the sensitivity of the film thickness and flow rate 
to inlet zone geometry was evaluated. The inlet type investigated was one 
where an inlet zone is induced with an expansion ring. Complete solutions 
have been obtained for several ring geometries by matching solutions in the 
three zones. The analysis predicts very strong sensitivity near the inlet 
that will be wear sensitive. The ring-induced inlet is predicted to provide a 
larger film and lower sensitivity to wear than does the premachined inlet 
currently incorporated in pumping ring seals. 

Experiments were performed on selected seal geometries for several seal mate- 
rials ranging in elastic modulus from 1.59 x 10^ MPa for a filled PTFE to 
4.58 X 10^ MPa for poly(amide-imide) . Two sets of tests were performed on 
each seal. One test set was conducted with the test seals flooded with oil on 
what would normally be the side exposed to the high-pressure gas. The second 
test set was performed under normal seal conditions with high-pressure gas on 
one side of the seal and ambient pressure oil on the opposite side. 

The flooded tests demonstrated the oil pumping capability of the experimental 
seals, a performance parameter readily calculated. The test conducted under 
normal seal operating conditions was used to demonstrate the gas sealing capa- 
bility of the expanding ring inlet PL seal design. I 

This report covers the PL seal analysis including a numerical design example, 
experimental results of several flooded and gas bearing tests, and a 
discussion of results. 
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2.0 EHD ANALYSIS OF HYDRODYNAMIC PUMPING RING 

2.1 The Mathematical Model 

A schematic of the Leningrader seal is shown in Figure 2-1. The seal is 
mounted on the shaft with an interference fit and is separated from the oil 
side by a backup spring acting on the secondary static seal as well as by the 
force generated by the high-pressure gas. A converging inlet region provides 
pumping action during the forward stroke. During the backstroke, the seal 
rubs and wipes the oil away. Any oil that leaks toward the high-pressure gas 
will be pumped back toward the oil reservoir. The pressure po at the inlet to 
the seal corresponds to a case where the seal is isolated from pg (by some 
secondary seal); if it is not, then pg = pg. While over the converging zone, 
the pressure will vary over most of the seal; past the initial point of 
contact, the pressures will be constant and equal to pg, the hoop stresses 
induced by shaft interference. 

The mathematical model of the seal is sketched in Figure 2-2. In part (a) the 
seal is portrayed as it would have looked shrunk onto the shaft without the 

I 

presence of an expansion ring; the dotted line A would be the seal's inner 
diameter with the shaft removed. When a ring is inserted at the outer end of 
the seal so as to stretch the ID of the seal, as shown in Figure 2-2(b), a 
flared opening is created at the end. Barring the presence of a lubricant, a 
reciprocating shaft would rub against the seal except over the ring-induced 
gap. Part (c) of Figure 2-2 portrays seal operation when a fluid film is pres- 
ent, and its hydrodynamic pressures separate shaft and seal during the forward 
stroke. The stresses and pressures corresponding to these three operational 
modes are shown in Figure 2-3. 

2.2 The Elasticity Equation 

Without the presence of a hydrodynamic film the configuration of the seal is 
determined approximately by the elasticity equation for an axisymmetric 
shell . 
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Fig. 2-2 


Mathematical Model of Leningrader Seal 


a) Interference — No Ring 
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c) Ring Present — Hydrodynamic Film 


Fig. 2-3 Stresses or Pressures in Mathematical Model 
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( 1 ) 


where the flexural rigidity D = Et ^ /12( l-V^ ) . With the nomenclature speci- 
fied in Figure 2-A, the boundary conditions are 


X = 0; 



t • • 

w = w = 0 


X "It » 



w = 0 


with the extra condition serving to determine L^, an unknown. The solution to 
equation (1) will be in the form of w = w (x), thus 6^ = 6^ (Lj-)* and the 
expression would have to be inverted to obtain in terms of the known input 

6 . 

r 


In dimensionless form the differential equation confined to the interval 0 < x 
< Lj- becomes 


d w 
dx^ 


+ w = 0 


0 < X < L, 


where 


w = - 



with <5 = 

o 


Et 


Po> 


( 2 ) 


The dimensionless coordinate x = -(x/H ) has reversed the direction of the 

e 

positive and negative x axes; and the normalizing factor is a diffusion 
length 



v^t 

[12 (1 - 


representing the extent over which an applied force or stress makes itself 
felt in the seal's shell-like structure. 


The boundary conditions in dimensionless form become 
X = 0; w = 1, w' = w" = 0 


(3a) 
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Fig. 2-4 Nomenclature for Seal Configurations 


r 


8 



X - ; w = (6^ + + 6^), w' = 0 


(3b) 


The above is a fourth order linear differential equation and its solution is 
of the form 


w 

where 


4 X.x 

E C. e ^ 

j=l ' 


X. 

3 


/2 



with the constant Cj determined by 


1 

X, 


Xie 




1 

X, 


X20 


^2^ 


1 

X, 


Xae 


X_L 
3 r 



Substituting in (4) w = 6^ we obtain 




4 X.(L /£ ) 

E C e ^ ® 

3 


(4) 


(5) 


which is the desired expression for Lf in terms of 6^. 

Figures 2-5 and 2-6 show the extent of lift-off in terms of both the ratio 

of the interferences (6/6 ) and their difference (6 - 6 ). It is seen that 

r s r s •_ 

most of the length of Lj. is formed at the very beginning. When = 26^, Lj. = 
2.75; when 6^ = 106^, Lj. = 3.9. Thus a fivefold increase in 6^ produces only 
a 42% rise in extent of Lj-. The shear on the other hand, except at the very 
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Ring Interference / Shaft Interference + 6^) / (S^ + S^) 


Fig. 2-5 


Relation Between Extent of Lift-off and Interferences 




Ring Location, 


Fig. 2-6 Extent of Lift-off as Function of Size of Ring 









beginning, rises about linearly with an increase in 6^, as shown in Figures 
2-7 and 2-8. At the beginning it behaves similarly to the lift-off, that is a 
small increase in 6^ produces a rapid rise in the shearing force at the point 
of contact x = 0. 

2.3 Elastohydrodynamic Solution 

2.3.1 The Governing Equations 

When there is a hydrodynamic fluid film h, its thickness, sketched in Figure 
2-9, is related to w via 

w = - (h + 6 ) (6) 

s 

where 6^ is a constant. The previously constant pressure Po in the gap is now 
replaced by the hydrodynamic pressure p(x). When these are introduced into 
equation (1), we obtain the elasticity equation in terms of h and p(x), namely 



The hydrodynamic pressures p(x) are obtainable from the one-dimensional 
Reynolds equation which, when integrated once, yields 



with ho a constant of integration to be determined. 

The boundary conditions for this EHD problem are 

X =-1. ;p = p,h=(5+5),h'=0 (9a) 

r g s r 

X ; P = P + ^ ^ ^o 

® R 

■*This term represents the hoop stresses in the seal which must be balanced by 
the hydrodynamic pressures. 
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^^8* 2—7 Shear Force at Contact Point as Function of Interferences 
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Fig. 2-9 Hydrodynamic Fluid Film Thickness 


In normalized form, the two differential equations and their boundary condi- 
tions become 


A 

dx^ 


+ h = p 


-B 

dx 


3 

(1 + h)-" 


— - *^3 + 

X = L ; - p = 1 +1 z ^ i » ^ 


(6 - 6 - h )/h = 0 

r s o o dx 


X lim p = 0, lim h = 0 


where 


e = 


6p u £ <s 

e s 



6p U Z 

e 



will be referred to as the EHD parameter. 


( 10 ) 


(11) 


(12a) 


(12b) 


Boundary conditions (12a) yield also the following useful relationships 



(1 + p) 


(13a) 


6 

r 


1 + h 

.1 + P 


(13b) 


6l ■ 


g 

(1 + p)^ 


(13c) 
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Other relationships between the normalized, quantities' are 


P - P (Lp 
1 + P (L^) 


h = 


1 + h 

1 + P ap 


p 



= p- p (ip = p + i + 





2.3.2 Solution of EHD Problem 


(13d) 


(13e) 


(13f) 


The origin of the coordinate system is here moved far upstream where h “ ho, 
as shown in Figure 2-10. Under these conditions h « 1 and so in the denomina- 
tor of equation (11), h can be ignored relative to 1. Dif f e.-entiating 
equation (10) and substituting for (dp/dx) from equation (11) we obtain 



dx 


Bh 


0 


(14) 


a linear equation in h (x). The solution to this differential equation is of 
the form 


h (x) 



m.x 


A.e ^ 
J 


where mj are the roots of the 5th order polynominal 
+ m + B = 0 


(15) 


Since h (x) 0 as x -►-<» we are allowed to retain only those roots that have 
positive real parts. It can be shown that equation (15) has only two such 
roots, a conjugate pair of complex numbers X and X. Combining these two roots' 
the solution then becomes 

h “ e Re. [e^^ (16) 
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where £ and (j> are real and X is one of the two roots of equation (15), namely 
the one having positive constants in both the real and imaginary parts. Thus 

h(0) = c cos <t> 


h^^(0)=eRe(X e*^) (n~th. derivative) 


In order to proceed with obtaining a solution, certain initial values are 
needed. We employ the following; 

• £ is assigned some small value, say 0.01 

• For any given B, X is determined from equation (15) 

• For any given value of <|> , values of h and its derivatives h'^^^ are 
determined for use as boundary conditions at x = 0 in equations (10) 
and (11). 

The above are used as starting points for a Runge-Kutta solution to equations 

A 

(10) and (11). The stepwise method is continued until (dh/dx) changes sign. 
We then 


• Obtain the corresponding value at Lj. by quadratic interpolation. 

• Carry out the Runge-Kutta integration over the range 0 - L^. 

• Obtain values of 5^ and B from equations (13b) and (13c), respective- 
ly, which in this form are functions of and B and thus must be 
inverted. However, it seems more efficacious to simply obtain direct 
computations of 5^ and Bj^ as functions of and B. 

2.3.3 Parametric Solutions 

The system under study here is well covered by the following ranges: 
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0.5 < < 50 
10 ^ < 6 < 10 ® 

The relations between these parameters and the variables <1> and are given in 
Figures 2-11 and 2-12. As seen, for large S, the 6^ = constant lines srs 
nearly parallel and we may thus approximate B by 

•.■Ct)"” 

where C is plotted in Figure 2-13. The deviation of B from B^ is shown in 
Figure 2-14. 

Single variable interpolation of the results shown in Figure 2-4 and bivariate 
interpolation of Figure 2-14 are used to compu te B . <l> is then computed from 8 
by a bivariate interpolation of the results in Figure 2-12. 

Figures 2-15 and 2-16 give the major items of interest in the EHD problem, 
namely the shapes and extents of the fluid film and its hydrodynamic pres- 
sures, for the case of Bj^ = 0.2. As seen, the parallel film portion is reached 
at a distance of 3 to 4 times 8-^. The pressures all peak, about the point where 
(dh/dx) -► 0 and are therefore located very near the point where contact would 
have occurred in the absence of a fluid film. As seen from Figure 2-17 the 
location of Pmax is very little affected by B and only mildly by 6^; a five- 
fold increase in 6^ produced a shift in of 20% for the 2 to 10 shift, and 
only 10% for the 10 to 50 shift in 6^. 

2.4 Solution with Seal Wear Included 
2.4.1 The Elasticity Equation 

It was postulated earlier, in Figure 2-2, that at the point of contact an 
infinite shear prevails. In practice this is unlikely as under such heavy 
loading the motion of the shaft would produce a high rate of wear at the 
initial point of contact and the shear concentration would be relieved. 
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This situation is portrayed graphically in Figure 2-18 where the shaded areas 
represent the amount of material removed by wear. Without wear, as postulated 
in previous sections, that amount of shaded seal substance would have had to 
be displaced above the shaft surface, thus causing the high shear concen- 
tration at the initial point of contact. Now, with that amount of interfer- 
ring matter removed, the stress is relieved so that, over the entire contact 
area, a uniform stress pg prevails. 

The differential equation now reads for the two separate regions bounded by 

X = 0 



p),-L <x<0 

o ’ r 

0 < X < “ 


(17a) 

(17b) 


The boundary conditions areJ 

(18a) 
(18b) 


X = ~Lj.» w = - 6^, w = 0 


X = 0; w, w , w , and w are continuous 


X = 00 


lim w = -6 


lim w 


II III 

w , w = 0 


(18c) 


Essentially the new condition introduced here is that the shear (w ) goes to 
zero at X = 0 along with the lower derivatives. 


In dimensionless form and with the reversed direction of the x axis the 
differential equation and the boundary conditions become 


.4 

d w 
dx^ 


+ w = 


0 , 


1 , 


X > 0 
X < 0 


(19a) 

(19b) 
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b) EHD Equivalent 


Fig. 2-18 The EHD Problem with Wear Included 
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(20a) 


— “ — ’ 

X = Lj.; w = 0 


_ _ _l —I I _• I t 

£ - 00 Lim w = 6, lim w , w , w = 0 


(20b) 


with aLL the derivatives continuous at x = 0. 


The solutions to equations (19a) and (19b) are 


4 

w = Z 

j=l 


t X . X 

C.e J 



Lj. > X > 0 


0 > X > - “> 


(21a) 


(21b) 


where Ot is an unknown constant. The four coefficients Cj and Ot can be 
obtained from use of the boundary conditions (20a) and (20b) yielding 



In the above, the zeros on the right-hand side of the matrix indicate not zero 
slopes but zero differences in slope at x = 0, in accordance with the require- 

I II III 

ment of continuity of w , w , and w at x = 0. 


After obtaining the values of Cj and Ot, the relation between Lj. and 6^ is 
obtained by first locating the value of x where w =0, finding the corre- 
sponding value of w which then supplies 5^ via 
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w = 


( 22 ) 


6 _ + 6 
^ o 

6+6 
s o 


1 + 6 ^ 


The shape of the separation gap due to seal wear in the vicinity of x =0 is 
portrayed in Figure 2-19. As seen, the extent of the lift-off gap had shrunk; 
and at the new point of contact the slope is no longer zero. The decrease in 
the separation length and the increase in slope at contact point with a 
rise in 6^ are given in Figures 2-20 and 2-21. 

2.4.2 The Hydrodynamic Solution 

The hydrodynamic solutions will here be approximated by representing the gap 
between the seal and shaft as equivalent to a composite axisymmetric bearing, 
as shown in Figure 2-18(b). The wedge-shaped portion of the bearing is given 
a slope equal to the slope of the seal at the initial point of contact with the 
shaft. The wedge can be made to extend to infinity, as with increasing h, the 
contribution to the pressures will decay rapidly. In the parallel portion of 
the bearing, the pressures must be constant and equal to pg, the hoop stresses 
in the seal. 


The details of solving this simple one-dimensional Reynolds equation with the 
proper boundary conditions are given in Reference 1. The solution for hg, 
which is the main quantity desired, is also given in Reference 1 as Equation 
A-15, namely 

. 3yU 

° Pg ■ Po^ ® 


Writing for the bearing slope 


6 = 


' 6+6 
s o I - 


w (0) 


we obtain for the value of the hydrodynamic gap 

3yU Z 


h = 


o (6 + 6 ) w' (0) (p + p - p ) 

o s g o 


(23) 
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Fig. 2-21 Slope at X = 0 for Worn Seal 







2.5 Numerical Computation 


Computations of displacements and pressures occurring in a Leningrader seal 
are performed for a set of typical conditions as given in Table 2-1. 

Figures 2-22 and 2-23 give the shapes and extents of the film thicknesses and 
pressure profiles for three levels of ring interference. The characteristics 
of the seal with the intermediate level of interference are summarized in 
Table 2-2. 

The effects of seal wear on hg as a function of ring interference 5^ is given 
in Figure 2-24. For the above case of = 0.5 mm, the 6.5 Um film thickness 
would upon wear increase to 7.8 pm. 

Although interference pressures, film thicknesses, and flow rates can be 
calculated based on the curves presented herein, a computer program PLGRAD has 
been provided for calculating these quantities for both worn and unworn seals. 
The program together with a sample input and output are provided the Appendix 
to this report. 
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Table 2-1 


Conditions for Numerical Example 


Outside Diameter of Seal, Dq, 

Inside Diameter of Seal, D{, 

Shaft Diameter, Dg, 

Viscosity, U 
Elastic Modulus, E, 

Speed, M, 

Poissons Ratio, v 
Stroke, S 

Shaft Radius, Rg = (Dq Di)/4 
Seal Thickness, t = (Dq - Di>/2 
Shaft Interference, * (Dg - Di)/2 
Shaft Pressure, pg = EtS/R^ 

^ ^ / A 

Diffusion Length = /^/[12(1 - v )] 


18.00 mm 
14.85 mm 

15.00 mm 
55 cps 

1.72375 GPa (0,25 x 10®psi) 
4000 rpm 
0.46 

34.00 mm 
8.2125 mm 

1.575 mm 
0.075 mm 

3,019 MPa (440 psi) 

2.051 mm 
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Table 2-2 


Results of Numerical Example 


Item 

SI 

Units 

English 

Units 

Ring Interference 

0.5 

mm 

20 

10"^ in 

Extent at Lift-Off 

7 

mm 

0.28 

in. 

Extent at Hydrodynamic Pressures 

10 

mm 

0.40 

in. 

Peak Pressure Location 

7.5 

mm 

0.30 

in. 

Peak Pressure 

7.7 

MPa 

1132 

psi 

Height of ho 

6.5 

)im 

0.28 

10'^ in 

Hoop Pressure, po 

3.2 

MPa 

470 

16 psi 
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3.0 EXPERIMENTAL PROGRAM 


The analytical treatment of ring-expanded pumping Leningrader (PL) seals 
provided direction to the task of designing actual functioning seals. An 
experimental testing program was then developed to provide data for the vali- 
dation of the analysis and to assess the ring-expanded PL seals performance as 
a working gas seal. 

The experimental test program was divided into two segments. The first, 
called the "flooded" series, tested PL seals in a total oil environment where 
oil at a slightly positive head was supplied to what is normally the gas side 
of the seal. The purpose of the flooded tests was to evaluate the oil pumping 
capabilities of the PL seal and to assist in the confirmation of the analyt- 
ical prediction for pumping rates and by inference to the film thickness pre- 
dictions. 

The second test series, referred to as the "gas" tests, evaluated the PL seals 
in an engine simulation environment with high-pressure gas on one side of the 
test seal and low pressure "crackcase" oil on the other. The results of these 
experiments, along with a description of the test apparatus, are discussed in 
the following sections. 

3.1 Test Apparatus 

The PL seal experimental program was carried out on reciprocating test appara- 
tus previously employed to evaluate hydrodynamic pumping rings (Ref. 1 and 2). 
The construction of the apparatus is modular and as such was readily modified 
to permit testing of pumping Leningrader seals. 

The test apparatus, designed and constructed to meet the requirements of this 

program, provides the following parameters: 

Rod Diameter: 19.0 mm 

Rod Stroke: 38.1 mm and 30.8 mm 

Rotational Speed: Variable over the range 10 < N < 60 Hz 
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Low-Pressure Oil Supply: 20W40 Oil at pressures to 0.41 MPa 
High-Pressure Gas Supply: He at pressures to 10.3 MPa 

The lower section of the apparatus contains the cranking mechanism illustrated 
in Figure 3-1. The crankcase assembly contains a crankshaft (1) that is 
supported on a pair of automotive-type sleeve bearings (2). An oil seal (3) 
is provided at the point where the crankshaft penetrates the crankcase for 
connection to a driver. The crankcase is split horizontally to allow for 
changing crankshafts. Removable counterweights (4) permit balancing of the 
crankshaft when changing either the reciprocating mass or the crank throw. 

The crankshaft assembly is bolted to a separate base that contains a varia- 
ble-speed electric motor. A gear tooth belt connects the drive motor to the 
crankshaft. Isolation mounts are attached underneath the base to provide 
vibration isolation. 

The test head components of the tester are directly connected to the crank- 
case. These components are identified on the assembly illustrated in Figure 
3-2. In this apparatus, a crosshead (1), driven by a connecting rod (2), 
applies reciprocating motion to the test rod (3). The crosshead is guided as 
it reciprocates by a crosshead sleeve (4). A lower housing (5), which serves 
as a foundation for the remainder of the apparatus, includes a manifolding 
system for providing oil to the crosshead sleeve. 

The upper end of the crosshead (1) contains a female taper, which is matched 
to the taper ground on the end of the test rod (3). Although the taper is 
sized for locking, a retaining clip (6) is also provided to assure that no 
relative axial motion exists between the crosshead and the test rod during 
tester operation. 

Mounted directly on the lower housing is a lower guide bearing housing (7), 
which is designed with a tight piloting arrangement to provide radial posi- 
tioning. The guide bearing (8) is an externally pressurized, pocketed design, 
which is supplied with the same oil as the crankcase, thereby eliminating the 
need for special seals. The oil supply to the guide bearing is adjustable to 
provide hydraulic shimming of the test rod for alignment purposes. Addi- 






PL Seal Housings 



Fig. 3-2 Test Head Assembly 
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tionaLly, the guide bearing functions to provide rod cooling which results 
from the oil exiting the bearing; this design eliminates the need for any 
auxiliary cooling systems. 

Tightly piloted and firmly attached to the lower guide bearing housing are two 
PL seal housings. A second bearing housing (9) containing a guide bearing (8) 
is tightly piloted to the upper PL seal housing to ensure proper bearing 
alignment. A simple cap (10) closes this upper bearing housing. 

The two housings and the PL seal installation is shown in Figure 3-3. For this 
configuration, two seals (1) oriented to pump oil away from each other are 
tested simultaneously. When fully assembled, the test rod (2) passes through 
the housing halves (3) which hold the seals (1) and act as a pressure vessel 
for the liquid or gas to be contained by the seals. Secondary seal seats (4) 
with integral 0-rings and the seal loading sleeves (5) provide the secondary 
seals for the PL test seals. Springs (6) provide an axial load to seat the 
seals. The thickness of shims (7) is adjusted to provide a range of prese- 
lectable axial loads on the seal. The separating plate (8) allows assembly of 
the housing halves with the seals installed. Outboard of each seal, annular 
grooves (9) in the housings supply oil for rod cooling and seal lubrication. 
Oil drains (10) empty the oil leakage gas separation cavity (11). PTFE 
inserts (12) maintain a close running gap for leakage containment without the 
possibilities of rod damage. 

For tests under flooded conditions, performed to corroborate analytical pre- 
dictions, oil at a positive head was supplied to the cavity between the two 
test seals in the area which is normally exposed to pressurized gas. For 
starved tests, high-pressure gas rather than oil was introduced into the same 
cavity so that the seals could be tested at conditions similar to those of rod 
seals in Stirling engines. 

Referring again to Figure 3-3, pressurized gas for gas sealing tests was 
introduced to the cavity at the center of the housing halves (3). During 
testing, any leakage gas (13) was flushed out of the separator cavity (11) by 
the introduction of a known amount of nitrogen gas. Any oil leakage was 
removed from the seal cavity through flushing ports (14). 
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3.2 Material Selection 


Material selection for the fabrication of PL test seals was based primarily on 
the need for candidate materials possessing low elastic modulus values. This 
requirement was necessary in order to permit reasonable manufacturing toler- 
ances on the seal bore diameter, the diameter which determines the interfer- 
ence fit between the seal and test rod. A second selection criterion was that 
the seal materials exhibit a resistance to cold flow and provide good dimen- 
sional stability. 

The final materials selected for PL test seals are listed in Table 3-1. The 
initial test material choice did not include the PTFE compound since that 
material had been extensively tested during related programs. Unsatisfactory 
test results on some of the primary materials, to be described in later report 
sections, led to inclusion of this material in the final test matrix. 

From the many available compounds, the materials indicated on Table 3-1 repre- 
sent a good variation in flexural modulus coupled with what was felt to be 
important secondary considerations including low friction coefficients, good 
dimensional stability (other than the PTFE compounds), and reasonable resist- 
ance to deterioration in an environment which could reach 125°C. 

3.3 Test Seal Geometry 

The design restriction placed on the PL seal geometry was to achieve minimum 
stress levels while maintaining both sufficient rod interference for proper 
pumping capacity and adequate expansion ring interference to promote a good 
inlet geometry, while meeting the maximum envelope restrictions imposed by the 
test vehicle. The result of these design requirements produced the PL seal 
envelope geometry illustrated by Figure 3-4. Table 3-2 lists the final seal 
bore dimensions and the actual interference values under which the seals were 
tested. 

The inlet geometry of the test seals is generated by the expansion of the seal 
body resulting from the expansion ring interference. Profilometer traces of 
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Table 3-1 


Pumping Leningrader Seal Material Selection 


Base Material 

Additives 

Manufacturer and 
Trade Name 

Elastic 

Modulus 

Modified Phenylene 
Oxide 

Unfilled 

General Electric Noryl SE-1 

2.48 X 

10^ 

MPa 

Polyetherimide 

Unfilled 

General Electric Ultem 1000 

3.00 X 

10^ 

MPa 

Po ly (amide- imide) 

Unfilled 

Amoco Chemicals Torlon 4203 

4.58 X 

10^ 

MPa 

PTFE 

Polyimide 

Powder 

Dixon Rulon J 

1.59 X 

10^ 

MPa 
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Measurements are in mm unless otherwise s 



Fig. 3-A Pumping 


“B" dia 


— 19.8 dia 
-28.6 dia 


eat for Ring Expander 
ore for Rod Interference 


jningrader Seal Geometry 
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Table 3-2 


Test Seal Interferences 

•Shaft Diameter: 19 mm nominal 

•Ring Expander Diameter: 20.6 mm nominal 


Seal 

Material 

Diametral* 
Shaft Interference 
10"° m 

Diametral* 
Ring Expander 
10~° m 

Modified Phenylene Oxide 

45.8 and 53.3 

206 and 211 

Polyetherimide 

50.8 and 48.3 

211 and 206 

Poly(amide-imide) 

50.8 and 48.3 

208 and 208 

PTFE 

140 and 137 

328 and 333 


*lst number is for upper test seal; 2nd number is for lower 
test seal 
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the deflected inlet shape .of the seals prior to testing are reproduced on 
Figure 3-5. The illustrated contours agree with predicted shapes. 

3.4 Test Method and Results 


The objective of the pumping Leningrader seal test program was to evaluate the 
performance of several seal designs over a range of operating conditions as a 
means of validating the advanced seal analysis. The fixed parameters under 
which the seals were tested include: 


•Oil Type: 

• Gas Type; 

• Inlet Oil Temp: 

• Reciprocating Stroke: 

• Rotational Frequency; 


20W40 

He 

49®C 

38.1 and 50.8 mm 
10, 35, 60 Hz 


Two test sequences were employed for experimentally verifying the PL seal 
analysis. Test Sequence 1, the "flooded" test in which oil at a positive 
pressure was supplied to the inlet of the test seals, was intended to provide 
the necessary data for correlation with analytically predicted flow charac- 
teristics. Test Sequence 2, the "gas" tests in which the test seals were 
evaluated in a more realistic environment with high-pressure gas, was intended 
as a verification of the ring-expanded PL seal design as a capable reciprocat- 
ing gas seal . 


Each of the two test sequences were evaluated under their own unique test 
conditions which included: 


• For "flooded" tests 

- oil inlet pressures of 0.21 and 0.41 MPa 

- oil discharge pressure of 0+ mm Hg 

- all seal materials tested at 50.8 mm stroke 

- one material tested at 25.4 mm stroke 

• For "gas" tests 

- gas inlet pressures of 10.3 and 5.2 MPa 

- oil discharge pressures of 0+ mm Hg and 0.069 Mpa 
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Relative Radial Dimension, mm 





Prior to actual seal testing, the level of preload exerted on the test seals 
by the preload springs was established. The proper seal load is achieved by 
adjustment of the shims controlling seal spring compression (Item 6 of Figure 
3-3). The shims were adjusted to provide a 365-N spring load for all the seal 
tests reported. This level of preload is sufficient to maintain secondary 
seal compression both for start-up conditions and for flooded tests. The high 
pressure exerted on the seals during gas sealing tests provided additional 
secondary seal loads. 

3.4.1 Flooded Tests 

The three primary materials were tested for pumping capability at the long 
stroke of 50.8 mm. The results of these tests are presented in Tables 3-3, 
3-4, and 3-5 for the modified phenylene oxide, the polyetherimide, and the 
poly(amide-imide) , respectively. Recorded on these tables is the test rod 
reciprocating frequency, oil inlet pressure and temperature, seal temper- 
atures, and the delivered flow (sum of two seals). 

Data points at the specified rod speeds reflect data taken over a long period 
of time. To assist in viewing these data for performance trends, the mean 
temperature levels and flow rates for each speed and inlet pressure combina- 
tion were examined. Although the variance of each mean temperature may appear 
excessive, these means, as presented in Table 3-6, do indicate trends not 
entirely evident when viewing the nonaveraged data. 

Just prior to the conclusion of the long-stroke tests on the modified pheny- 
lene oxide test seals, erratic flow behavior was noted. Examination of both 
test seals revealed longitudinal cracks had occurred in the thin wall sections 
extending from the expansion ring to the start of the trapezoidal cross 
section. Since the seals had already been designed to provide minimum stress 
levels for the available envelope restrictions, it was concluded that testing 
of the modified phenylene oxide in its present configuration should be discon- 
tinued. 

A replacement material of PTFE with a polymer powder additive was selected for 
gas seal testing so that a three-material test matrix could be maintained. 


53 



Table 3-3 


Flooded Test Results - Modified Phenylene Oxide 
Stroke 50.8 mm. Oil = 20W40 (Mobil Delvac) 


Speed 

cpm 

Inlet 
P MPa 

■H 





10 

0.21 

30 

48 

49 

420 


10 

0.21 

29 

48 

48 

429 

(426) 

10 

0.21 

29 

47 

47 

429 


10 

0.41 

30 

53 

54 

411 


10 

0.41 

30 

51 

52 

438 

(432) 

10 

0.41 

31 

49 

49 

447 


35 

0.21 

26 

48 

46 

742 


35 

0.21 

28 

54 

54 

742 

(739) 

35 

0.21 

31 

59 

53 

733 


35 

0.41 

29 . 

114 

52 

769 


35 

0.41 

30 

59 

59 

733 

(748) 

35 

0.51 

31 

60 

60 

742 


60 

0.21 

29 

52 

52 

1317 


60 

0.21 

34 

64 

63 

1218 

(1251) 

60 

0.21 

37 

69 

68 

1218 


60 

0.41 

37 

75 

71 

1152 


60 

0.41 

38 

75 

72 

1053 

(1075) 

60 

0.41 

38 

77 

73 

1020 



(1) Upper Seal 

(2) Lower Seal 
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Table 3-4 


Flooded Test Results - Polyetherimide 
Stroke 50.8 mm, Oil = 20W40 (Mobil Delvac) 


Speed 

Inlet 

Inlet 

Temp. 

Temp. 

Flow 


Hz 

P MPa 

T °C 

“C(D 

°C(2) 

cc/min , 

(Average) 

10 

0.21 

28 

47 

45 

250 


10 

0.21 

28 

48 

48 

277 

(271) 

10 

0.21 

28 

48 

47 

286 


10 

0.41 

29 

58 

56 

260 


10 

0.41 

29 

53 

55 

259 

(262) 

10 

0.41 

29 

53 

52 

268 


10 

0.41 

29 

52 

50 

259 


35 

0.21 

28 

52 

52 

465 


35 

0.21 

28 

56 

53 

. 554 

(521) 

35 

0.21 

28 

60 

58 

545 


35 

0.21 

26 

43 

43 

554 


35 

0,21 

27 

49 

51 

572 

(560) 

35 

0.21 

28 

58 

60 

554 

35 

0.41 

28 

61 

63 

554 


35 

0.41 

28 

64 

67 

554 

(554) 

60 

0.21 

28 

60 

64 

751 


60 

0.21 

29 

76 

77 

662 

(694) 

60 

0.21 

29 

79 

80 

670 


60 

0.41 

29 

83 

84 

670 


60 

0.41 

29 

87 

87 

670 

(672) 

60 

0,41 

29 

89 

88 

679 


60 

0.41 

30 

89 

86 

670 



(1) Upper Seal 

(2) Lower Seal 
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Table 3-5 


Flooded Test Results - Poly(amide-imide) 
Stroke 50.3 mn, Oil = 20M40 (Mobile Delyac) 


Speed 

Hz 

Inlet 
P MPa 

Inlet 

T.°C 

Temp. 

“C(D 

Temp 

°C(2) 

Flow 

cc/min 


0.21 

28 

48 

49 

322 

10 

0.21 

29 

49 

50 

331 

10 

0.21 

29 

50 

51 

322 


0.41 

30 

64 

66 

304 

10 

0.41 

30 

60 

62 

313 

10 

0.41 

30 

45 

58 

304 

35 

0.21 

26 

47 

42 

626 

35 

0.21 

27 

60 

81 

608 

35 

0.21 

28 

66 

67 

608 

35 

0.41 

28 

69 

68 

617 

35 

0.41 

29 

72 

70 

608 

35 

0.41 

29 

73 

72 

617 

60 

0.21 

29 

68 

66 

769 

60 

0.21 

30 

79 

76 

706 

60 

0.21 

30 

85 

82 

688 

60 

0.42 

31 

86 

83 

679 

60 

0.42 

32 

89 

87 

670 

60 

0.42 

32 

90 

89 

688 


(1) Upper Seal 

(2) Lower Seal 
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Table 3-6 

Mean Temperature Levels and Flow Rates 


Seal 

N 

P 

Tmean 

AT 

Flow 

Material 

Hz 

MPa 

“C 

°C 

cc/min 

Modified 

10 

0.21 

48/48 

18/18 

426 

Phenylene 

35 

0.21 

54/51 

26/23 

739 

Oxide 

60 

0.21 

62/61 

28/28 

1251 

(Noryl) 

10 

0.41 

51/52 

21/22 

432 


35 

0.41 

58/57 

28/27 

748 


60 

0.41 

76/72 

39/36 

1075 


Polyetherimide 10 

0.21 

47/47 

19/18 

271 

^ 35 

0.21 

50/51 

23/24 

560 

(Ultem) gQ 

0.21 

72/74 

43/45 

694 

10 

0.41 

54/53 

24/24 

261 

35 

0.41 

63/65 

35/37 

554 

60 

0.41 

87/86 

58/57 

672 


PolyCamide- 

10 

0.21 

imide) 

35 

0.21 

(Torlon) 

60 

0.21 


10 

0.41 


35 

0.41 


60 

0.41 


49/50 

20/22 

325 

57/56 

30/28 

614 

78/74 

48/45 

721 

60/62 

30/32 

307 

71/70 

43/42 

614 

88/86 

57/54 

679 






A short-stroke (25.4 mm) test was run using the polyetherimide seals. At this 
stroke, which was less than the overall length of the seals themselves, oil 
pumping occurred only at maximum conditions of speed and inlet oil pressure 
and then only as miniscule amounts. After several unsuccessful attempts to 
improve the flow rates, the short-stroke tests were terminated. 

3.4.2 Gas Sealing Tests 

At the conclusion of the flooded tests, the apparatus was refitted for gas 
sealing tests. At the same time the new set of PTFE seals was inspected for 
rod and expansion ring interference. Because of the known tendency for PTFE 
to exhibit permanent dimensional changes when under stress, the nominal inter- 
ference levels were increased above that of the other seals. 

The rod interferences for the PTFE seals were measured at 140 and 136 x 10 ® m 
for seals 1 and 2, respectively. The expansion ring interferences for the 
same seals were 328 and 332 x 10”^ m. 

The flooded tests had confirmed the high forward stroke pumping capability of 
the ring-expanded seals. It was important for the gas test however to assure 
that some oil for seal lubrication would be brought into the seal-rod inter- 
face on the return stroke. The oil entering the gas side of the seal would 
then be removed by the strong pumping action of the ring-expanded seal's geom- 
etry. 

To provide a weak but otherwise effective pumping action on the return stroke, 
each test seal was provided with a small inlet chamfer as illustrated by 
Figure 3-6. 

Actual chamfer dimensions were confirmed by profilometer trace measurements 
and are shown in Table 3-7. 
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Forward 

Stroke 

Inlet 

(Gas Side) 



Reverse 

Stroke 

Inlet 

(Oil Side) 


Oil Inlet 
Chamfer 


Fig. 3-6 Test Seal with Small Inlet Chamfer 
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Table 3-7 


Return Stroke 

Inlet Chamfer 

Geometry 

Seal Material 

Chamfer < 

Chamfer Extent 

Polyetherimide 

0“ - 31' 

0.145 

Poly ( ami de-imide) 

0° - 35’ 

0.148 

PTFE 

1° - 02* 

0.150 


Gas sealing tests were conducted for various combinations of gas pressure, oil 
inlet pressure, and reciprocating frequency at one stroke. For comparison 
purposes, one set of test conditions was run at a shorter stroke. The follow- 
ing list indicates the values of the parameters actually tested: 

Gas Pressure: 5.2 and 10.3 MPa 

Oil Pressure: 0+ and 0.07 MPa 

Reciprocating Frequency: 16.7, 30.0, 60.0 Hz. 

Stroke: 50.8 iran primary, 38.1 mm secondary 

The three seal materials were tested at all combinations of pressure (both gas 
and oil) and frequency. For each combination of test conditions, the seals 
were run until temperature stabilization occurred, at which time any gas leak- 
age was measured and the seal temperatures recorded. To end each test cycle, 
a raid-speed and lower gas pressure combination was run for several hours. 

At the completion of each test point the tester was disassembled and oil 
wetted parts in the gas side of the seal were examined and any residual oil 
weighed . 

The results for the poly(amide-imide ) seal tests are given on Tables 3-8 and 
3-9, the polyetherimide seal tests on Tables 3-10 and 3-11, and the PTFE seal 
tests on Tables 3-12, 3-13, and 3-14. 
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Table 3-8 


Poly(amide-imide) Gas Seal Test Results 

(Oil Leakage - Average for Run 0.51 g/hr; Total Run Time 3:30) 
(Stroke 50.8 mm; Oil Side Pressure 0+) 


Upper Lower 


Gas Side 
Pressure 
MPa 

Reciprocating Upper Seal 

Frequency Leakage Rate 

Hz cc/min 

Lower Seal 
Leakage Rabe'> 
cc/min 

Oil Inlet 
Temp 
°C 

Seal 

Temp 

“C 

Seal 

Temp 

°C 

5.2 

16.7 

0.15 

0.12 

36.8 

0) 

> 

59.6 

5.2 

30.0 

3.00 

0.45 

37.7 

•H 

U 

CO 

u 

58.5 

5.2 

60.0 

4.20 

4.80 

39.8 

CD 

a 

o 

62.8 

10.3 

16.7 

Off Scale 

0.36 

42.3 

M 

CD 

fH 

62.8 

10.3 

30,0 

Off Scale 

1.50 

42.7 

Cl 

3 

O 

O 

65.1 

10.3 

60.0 

1.65 

4.80 

42.9 

o 

a 

u 

(U 

66.3 

5.2 

30.0 

0.51/4.50 

1.35/4.80 

47.6 

JC 

H 

58.7 
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Table 3-9 


Poly(amide-imide) Gas Seal Test Results 


(oil Leakage - Average for Run 0.3 g/hr; Total Run Time 4:00) 
(Stroke 50.8 mm; Oil Side Pressure 0.07 MPa) 


Gas Side 
Pressure 
MPa 

Reciprocating 

Frequency 

Hz 

Upper Seal 
Leakage Rate 
cc/min 

Lower Seal 
Leakage Rate 
cc/min 

Oil Inlet 
Temp 
“C 

Upper 

Seal 

Temp 

°C 

Lower 

Seal 

Temp 

°C 

5.2 

16.7 

2.40 

0.26 

27.0 

0) 

> 

59.6 

5.2 

30.0 

0.30 

1.05 

30.6 

38.9 

5.2 

60.0 

Off Scale 

3.30 

31.6 

cO 

<U 

P- 

O 

44.3 

10.3 

16.7 

Off Scale 

6.60 

38.8 

c 

M 

46.9 

10.3 

30.0 

Off Scale 

6.30 

37.1 

0) 

Q- 

0 

48.0 

10.3 

60.0 

Off Scale 

6.00 

33.4 

0 

o 

0 

B 

49.7 

5.2 

30,0 

0.012/ 
Off Scale 

0.03/0.51 

34.2 

u 

<3) 

H 

47.8 
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Table 3-10 


Polyetherimide Gas Seal Test Results 



(Oil Leakage - Average for Run; 

(Stroke 50.8 mm; Oil 

; Trace; Total 
Side Pressure 

Run Time 
0+) . 

7:00) 


Gas Side 
Pressure 
MPa 

Reciprocating 

Frequency 

Hz 

Upper Seal 
Leakage Rate 
cc/min 

Lower Seal 
Leakage Rate 
cc/min 

Oil Inlet 
Temp 
“C 

Upper 

Seal 

Temp 

°C 

Lowe 

Seal 

Temp 

°C 

5.2 

16.7 

Off Scale 

13.8 

47.1 

37.1 

53.4 

5.2 

30.0 

0.69 

0.63 

46.7 

41.4 

56.6 

5.2 

60.0 

Off Scale 

5.10 

46.5 

47.0 

70.0 

10.3 

16.7 

19.5 

0.63 

47.2 

40.3 

61.1 

10.3 

30.0 

2.55 

2.55 

47.9 

41.2 

71.0 

10.3 

60.0 

3.00 

7.20 

48.8 

62.0 

74.2 

5.2 

30.0 

0.30/18.0 

1.32/2.70 

47.8 

57.5 

68.7 
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Table 3-11 


Polyetherimide Gas Seal Test Results 



(Oil Leakage - 
(Stroke 

Average for Run 2 g/hr; Total 
50.8 mm; Oil Side Pressure 0. 

Run Time 
,07 MPa) 

2:18) 


Gas Side 
Pressure 
MPa 

Reciprocating 

Frequency 

Hz 

Upper Seal 
Leakage Rate 
cc/min 

Lower Seal 
Leakage Rate 
cc/min 

Oil Inlet 
Temp 
“C 

Upper 

Seal 

Temp 

°C 

Lower 

Seal 

Temp 

“c 

5 .2 

16.7 

0.42 

1.05 

26.9 

38.7 

50.2 

5.2 

30.0 

0.90 

2.70 

28.2 

41.0 

54.0 

5,2 

60.0 

0.45 

1.50 

35.0 

51.4 

62.0 

10.3 

16.7 

10.5 

1.08 

36.2 

59.0 

65.1 

10.3 

30.0 

10.5 

1.20 

38.0 

64.8 

72.9 

10.3 

60.0 

Off Scale 

3.00 

39.4 

69.6 

80.8 

5.2 

30.0 

0.09/12.6 

0.27/0.81 

47.2 

58.8 

67.1 
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Table 3-12 


PTFE Gas Seal Test Results 

(Oil Leakage - Average for Run 0.3A g/hr; Total Run Time 5:00) 
(Stroke 50.8 mm; Oil Side Pressure 0+) 


Gas Side 
Pressure 
MPa 

Reciprocating 

Frequency 

Hz 

Upper Seal 
Leakage Rate 
cc/min 

Lower Seal 
Leakage Rate 
cc/min 

Oil Inlet 
Temp 
“C 

Upper 

Seal 

Temp 

“c 

Lower 

Seal 

Temp 

°C 

5.2 

16.7 

0.33 

0.27 

35.9 

30.8 

39.7 

5.2 

30.0 

2.70 

0.45 

37.3 

34.6 

46.2 

5.2 

60.0 

2.40 

3.30 

39.1 

43.1 

59.6 

10.3 

16.7 

0.45 

0.72 

42.3 

42.7 

56.3 

10.3 

30.0 

2.40 

0.54 

41.5 

42.3 

58.4 

10.3 

60.0 

2.40 

2.70 

39.7 

41.8 

66.0 

5.2 

30.0 

0.12/3.00 

0.42/2.70 

48.9 

46.3 

62.5 
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Table 3-13 



(Oil Leakage - 
(Stroke 

PTFE Gas Seal Test Results 

Average for Run, Trace; Total 
50.8 mm; Oil Side Pressure 0. 

Run Time 
,07 MPa) 

5:00) 


Qas Side 
Pressure 
MPa 

Reciprocating 

Frequency 

Hz 

Upper Seal 
Leakage Rate 
cc/min 

Lower Seal 
Leakage Rate 
cc/min 

Oil Inlet 
Temp 
“C 

Upper 

Seal 

Temp 

“c 

Lower, 

Seal 

Temp 

°C 

5.2 

16.7 

0.36 

0.10 

41.3 

33.6 

38.4 

5.2 

30.0 

0.22 

0.60 

42.0 

37.0 

40.8 

5.2 

60.0 

3.30 

3.60 

43.8 


48.2 

10.3 

16.7 

Off Scale 

0.09 

44.1 1 


47.8 

10.3 

30.0 

Off Scale 

2.10 

43.5 


49.2 

10.3 

60.0 

0.66 

5.70 

42.8 


53.3 

5.2 

30.0 

Off Scale 

0.04/6.30 

48.6 

55.2 

55.7 
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Table 3-14 


PTFE Gas Seal Test Results 



(Oil Leakage - 
(Stroke 

Average for Run, Trace; Total 
38.1 mm; Oil Side Pressure 0. 

Run Time 
07 MPa) 

3:30) 


Gas Side 
Pressure 
MPa 

Reciprocating 

Frequency 

Hz 

Upper Seal 
Leakage Rate 
cc/min 

Lower Seal 
Leakage Rate 
cc/min 

Oil Inlet 
Temp 

“c 

Upper 

Seal 

Temp 

“C 

Lower 

Seal 

Temp 

“c 

5.2 

16.7 

0.36 

0.10 

30.0 

30.4 

29.1 

5i2 

30.0 

0.06 

0.33 

33.0 

37.0 

36.8 

5.2 

60.0 

0.05 

0.69 

37.8 

46.0 

50.5 

10.3 

16.7 

0.07 

0.54 

41.0 

46. L 

44.2 

10.3 

30.0 

0.02 

0.13 

41.5 

50.0 

50.0 

10.3 

60.0 

1.68 

0.69 

46.3 

59.7 

60.2 

5.2 

30.0 

0.04/17.1 

'vO/1.56 

51.6 

56.6 

57.0 

5.2* 

16.7 

0.09/5.1 

2.40/8.70 

50.3 

53.6 

46.6 


*Data taken following day 
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4.0 DISCUSSION OF RESULTS 


4.1 Flooded Tests 

The purpose of conducting the flooded PL seal tests was to generate actual 
flow data which could then be used for correlation with analytically predicted 
performance. The flooded tests confirmed the excellent pumping capability of 
the expanded-ring inlet geometry of the new PL seal design; an important 
attribute if the seal is to successfully purge any oil which may enter the gas 
side of the seal . 

To illustrate the degree of correlation between the test data and analytically 
predicted flow values, each data item listed in Table 3-6, which represents 
the sum of flow rates for double-acting seals, is presented for comparison 
purposes with the flow rates calculated using PLGRAD (See Appendix) in Table 
4-1, 

In several areas, agreement between prediction and test is quite good, with 
major deviations occurring with underpredictions at 10 Hz and for the highest 
elastic modulus seal. 

The test data is consistent with expected results in that supply pressure has 
little or no effect on flow rate. There is evidence, as seen in Table 3-6, 
that flow rates even decrease slightly with increases in supply pressure. 
This decrease, however, may be more a result of changes in seal loading rather 
than from other, more obscure reasons. 

In addition, there seems to be only a modest effect on flow rates from vari- 
ations in seal operating temperatures. At each test frequency, data were 
collected over an extended period of operational time in order to monitor seal 
temperature stabilization. The effect on flow of an increasing seal temper- 
ature during an otherwise stable test condition is not zero but is definitely 
a minor consideration, at least over the range of temperatures experienced. 
This minimal effect is evident by the narrow scatter band of pumped flow rates 
at each fixed speed and inlet pressure. 
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Table 4-1 


Experimental and Analytical Flow Rate Comparison 
Flooded Tests 


Seal 

Inlet P 

M 

Tmean 

Flow Rate 

cc/min 

Material 

MPa 

Hz 

“C 

Experiment 

Analysis 

Modified 


10 

48/48 

426 

64.2 

Phenyl ene Oxide 

0.21 

35 

54/51 

739 

656.0 



60 

62/61 

1251 

1356.0 

(Noryl) 

— 

10 

51/52 

432 

57.1 


0.41 

35 

58/57 

748 

549.0 



60 

76/72 

1075 

942.0 


Polyetherimide 


10 

47/47 

271 

49.8 


0.21 

35 

50/51 

560 

489.0 

(Ultem) 

— 

60 

72/74 

694 

858.0 



10 

54/53 

261 

44.4 


0.41 

35 

63/65 

554 

382.0 



60 

87/86 

672 

743.0 


Poly(amide-imide ) 


10 

49/50 

325 

20.9 

0.21 

35 

57/56 

614 

200.0 

(Torlon) 

— 

60 

78/74 

721 

364.0 



30 

60/62 

307 

13.6 


0.41 

35 

71/70 

614 

132.0 



60 

88/86 

679 

245.0 
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There are two other parameters affecting pumped flow rates. One is the obvi- 
ous speed dependent effect, and the other, a more subtle effect resulting from 
material properties. 

As speeds were increased above the first test speed of 10 Hz, a rapid rise in 
pumped flow rate was experienced. As the maximum test speed of 60 Hz was 
approached, however, the rate of oil pumping seemed to be approaching a maxi- 
mum. The material effect is demonstrated by the much higher pumping flow rate 
exhibited by the Phenylene Oxide seals with lower flow rates produced by the 
higher modulus material. 

The speed effect is affected by two processes acting in concert. As recipro- 
cating rod velocities are increased, more fluid per unit time can enter the 
seal inlet to generate higher film pressures and a larger film thickness, 
thereby producing an increased flow rate. At the same time, however, at the 
transition between the forward and rearward rod velocity directions, the elas- 
ticity of the seal must collapse the fluid film in order to shut off any 
reverse flow. At very low speeds the oil volume brought into the inlet is low 
and the time available to collapse the fluid film by squeeze-film effects is 
high, thereby providing a low flow rate with minimal backflow. At moderate 
speeds the volume of oil brought into the seal is increased but some backflow 
exists since the fluid squeeze film collapse time is now a larger percentage 
of the stroke period. At higher speeds, although the volume of oil introduced 
into the seal is large, the squeeze film collapse time is now a significant 
portion of the cycle time, and backflow flow rates can be a significant 
portion of incoming rates; the net result being a leveling off of the pumped 
flow rates. 

The effect of PL seal material elastic modulus on the rate at which seals pump 
is attributed to restrictions to film thickness development caused by higher 
elastic modulii. Lower values of film thickness produced by higher material 
stiffness properties would inhibit the development of high flow rates. 

Two test results differ significantly from analytical predictions; one is the 
high flow rates produced by all seals at a speed of 10 Hz, the second is the 
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higher than predicted flows produced by the poly(amide”imide) seals. Both 
results may be explained as follows. 

At 10 Hz, where sufficient time at the end-of-stroke direction would permit 
fluid film collapse, the static friction force developed between the seal and 
the test rod could be sufficient to diminish the effectiveness of the PL 
seal's outer diameter static seal (Figure 2-1). This loss of sealing capacity 
could allow oil to by-pass the PL seal and bias its actual performance, there- 
by producing the results seen. Since the preload spring for the test seals 
has a very high spring rate (1.74 x 10^ N/M), very small seal motions would 
quickly overcome any static friction force development to limit any lift-off; 
this, combined with a long drain path for the leakage flow, would provide 
sufficient flow resistance to prevent an extremely large amount of leakage to 
occur. 

The same arguments can be applied to the test results for the 
poly(amide-imide) test seals which show a much higher than anticipated flow 
rate. For this material, because of its higher modulus, sufficiently large 
friction forces could be developed at all speeds to provide a flow bias for 
all the test points. 

4.2 Gas Sealing Tests 

The results of the gas sealing tests for each of the three seal materials show 
wide variations in gas leakage at essentially similar test conditions. This 
variability even extends to the long-term tests where independent variables 
were maintained constant. Greater gas leakage variability does however seem 
to be associated with the seal materials with high elastic modulus values 
although even the PTFE seals suffer from this apparent phenomena. 

The accuracy of the gas leakage values depends on exactly how the gas (in this 
case He) which actually leaks past the seal mixes with the purged gas. If the 
purged gas is not well mixed with the leakage gas at the time of sampling, the 
helium leak detecting system may record higher proportional levels of helium 
in the purged gas mixture than are actually present. This problem seems to 
have occurred with some frequency since it was possible under certain circum- 


71 



stances to obtain large differences in detector readings over short time peri- 
ods at supposedly steady-state operation. 

One conclusion which can be readily drawn from the gas tests is that the wide 
variability in the data taken for each seal requires that large test samples 
of many seals must be evaluated before a statistically sound judgment of 
performance can be made. 

The evidence is strong that low modulus PL seal materials such as PTFE out- 
perform higher modulus materials such as poly(amide-imide) with respect to 
both oil and gas leakage. This conclusion is based, however, on single sample 
test results and may not hold up under the scrutiny of multiple sample test- 
ing. 

One good argument for the use of lower modulus materials for the construction 
of PL seals is the freedom this material type permits in establishing seal 
dimensions. Larger seal-to-piston-rod interferences permitted by low elastic 
modulus materials means than larger tolerance can be applied to the seal bores 
without concern for substantially lowering the interference contact pressure 
that would result from the same tolerance applied to stiffer materials. 

At the conclusion of the polyethermide test sequence, a visual examination of 
the test seals revealed that both the upper and lower seals had experienced 
serious structural failures. These failures manifested themselves as circum- 
ferential cracks in the trapezoidal sections of the seals; cracks which 
extended from the secondary seal surface axial through a portion of the heavy 
trapezoidal seal section. There is no evidence that these cracks, most likely 
the result of elevated shear stresses induced by the high sealed gas pres- 
sures, projected through the body of the seal to permit leakage of either gas 
or oil. The fact that the cracks appeared after a relatively short testing 
period (9 hr 18 min) does not provide confidence for the materials long-term 
durability. 

Of the two remaining materials which did not show distress after testing, the 
poly(amide-imide) with its very high elastic modulus still presents problems 
for the PL seal designer. The high modulus of this material requires that its 
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rod interference be kept small in order to minimize starting friction problems 
and to alleviate assembly difficulties which would result from forces gener- 
ated by large interference fits. 

The combination of the low interference fit and high elastic modulus of the 
poly(amide-imide) material also means less conformity of the material to the 
small irregularities which occur during manufacture. Relative lack of 
conformability may explain why this particular material showed high leakage of 
both gas and oil during the sealing tests. 
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5.0 DISCUSSION AND CONCLUSIONS 


The analysis presented in this report is capable of predicting seal perform- 
ance and sensitivity of the seal performance to wear. The ring-induced inlet 
appears to have low sensitivity to the wearing-in process and provides excel- 
lent pumping capacity and large films. 

Experimental work has demonstrated the ring-induced inlet PL seal pumps oil at 
high flow rates and has the ability to provide excellent sealing capability. 
Limited seal tests, however, did not permit a rigorous statistical evaluation 
of seal performance. 

The experimental performance and analytical performance predictions, along 
with design requirements for producing practical PL seals, result in the 
following conclusions: 

1. Low modulus materials such as PTFB, although possessing limited high 
temperature capability and cold flow resistance, produced the best 
PL seal performance. This performance advantage is attributed to the 
inability of high modulus materials to conform to slight machining 
irregularities, resulting in excessive gas leakage. The high elas- 
tic modulus materials also require very small rod interferences if 
the seal is to be assembled. This means that seals must be manufac- 
tured to more precise tolerances, always difficult in plastic. 

2. The ring-expanded inlet PL seal was shown to be an excellent oil 
pump. Predicted pumping rates, based on film thickness calcu- 
lations, are in reasonable agreement. 

3. The prediction of lower-than-experienced pumping rates is attributed 
to a lower prediction of velocity effects and the absence of 
squeeze-film behavior in the analysis. The squeeze-film analysis 
was not included in the present work since it is felt that a starved 
PL seal acting in a gas, oil-spray environment would exhibit little 
squeeze-film behavior. 
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since PL seals are not contemplated to be used as oil pumps, the addition of 
squeeze-film behavior to the present analysis should not have a high priority. 
It can, however, have a beneficial effect on seal lubrication which- should be 
considered. Of primary importance to seal development is the search for low 
modulus seal materials that do not exhibit cold flow from which better PL 
seals could be manufactured. With a wider selection of available materials, a 
statistically sound test program with a sufficent sample size should be initi- 
ated. 
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APPENDIX 


COMPUTER PROGRAM PLGRAO 

Computer program PLGRAD, which has been used to generate performance data used 
in this report, is listed herein. The FORTRAN listing of program PLGRAD 
FORTRAN is preceded by an input .description, an output description, and sample 
inputs and outputs. The sample input file and corresponding output file are 
for an unworn and worn phenylene oxide seal. 

A.1 Input Description 

The program uses namelist input. Definitions of the namelist quantities are 
given below. 


Parameter 
X No. X 

Namelist 
X Name x 

Definition 

1 

DIMM 

Seal ID (mm) 

2 

DSMM 

Shaft Diameter (mm) 

3 

DOMM 

Seal OD (mm) 

4 

STMM 

Stroke (mm) 

5 

VI CP 

Viscosity (Cp) 

6 

EPSI 

Elastic Modulus (Ib/in.^) 

7 

RPM 

Speed (rpm) 

8 

RIMM 

Radial Ring Interference (mm) 

9 

HPMP 

External-Inlet Pressure (MPa) 

10 

POIS 

Poisson’s Ratio 


IWORN 

NP 

DX 

0, Unworn or 1, Worn 
Parameter No. to be Varied 
Parameter Increments 


N Number of Points 

The program performs N calculations starting with the input value of parameter 
NP, then increments by DX, N-1 times. 

A.2 Output Description 

The outputs are in the form of a four-column table. The first column contains 
the parameters selected by the input value NP. The remaining columns are: 

HMIC = Film Thickness (10 ^ m) 

PMPA = Interface Pressure (MPa) \ 

CCPM = Flow Rate for Double-Acting Seal (cc/min) 

A.3 Sample Input 

& INPUTS 

DIMM=18.850, DSMM=19.00, DOMM=22.00, STMM=50.80, 

VICP=55., EPSI=0.38E6, RPM=1000., RIMM=.25, 

HPMP=0., POIS=0.46, IWORN=0, 

NP=7, DX=100., N=10, 

&END 

&INPUTS IWORN=l,&END 
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A.4 Sample Output 


RPM 

HMIC 

PMPA 

CCPM 

l.OOOOOD+03 

5.21409D+00 

2.96753D+00 

9.75272D+00 

l.lOOOOD+03 

5.57180D+00 

2.96753D+00 

1.14640D+01 

1.20000D+03 

5.91862D+00 

2.96753D+00 

1.32846D+01 

1.30000D+03 

6.25562D+00 

2.96753D+00 

1.52111D+01 

1.40000D+03 

6.58372D+00 

2.96753D+00 

1.72404D+01 

1.50000D+03 

6.90366D+00 

2.96753D+00 

1.93695D+01 

1.60000D+03 

7.21609D+00 

2.96753D+00 

2.15958D+01 

1.70000D+03 

7.52160D+00 

2.96753D+00 

2.39170D+01 

1.80000D+03 

7.82067D+00 

2.96753D+00 

2.63308D+01 

1.90000D+03 

8.11374D+00 

2.96753D+00 

2.88351D+01 


RPM 

HMIC 

PMPA 

CCPM 

l.OOOOOD+03 

5.96042D+00 

2.96753D+00 

1.11487D+01 

l.lOOOOD+03 

6.55646D+00 

2.96753D+00 

1.34899D+01 

1.20000D+03 

7.15251D+00 

2.96753D+00 

1.60541D+01 

1.30000D+03 

7.74855D+00 

2.96753D+00 

1.88413D+01 

1.40000D+03 

8.34459D+00 

2.96753D+00 

2.18514D+01 

1.50000D+03 

8.94063D+00 

2.96753D+00 

2.50846D+01 

1.60000D+03 

9.53667D+00 

2.96753D+00 

2.85407D+01 

1.70000D+03 

1.01327D+01 

2.96753D+00 

3.22197D+01 

1.80000D+03 

1.07288D+01 

2.96753D+00 

3.61218D+01 

1.90000D+03 

1.13248D+01 

2.96753D+00 

4.02468D+01 


A.5 PLGRAD FORTRAN Listing 


C$ H2 

C THIS IS THE DIMENSIONAL INLET PROGRAM. 

COMMON/BINP/DIMM,DSMM,DOMM,STMM,VICP,EPSI,RPM,SDEG,HPMP,POIS,DMMM 
C HERE HPMP DENOTES DIFFERENCE BETWEEN EXTERNAL AND INLET PRESSURES. 
DIMENSION A(ll) 

EQUIVALENCE (A(1),DIMM) 

NAMELI ST/ INPUTS / DIMM , DSMM , DOMM , STMM , 

+VI CP , EPS I , RPM , RIMM , HPMP , PO I S , DMMM , I WORN , NP , N , DX , I PAR 
INTEGER STR(11)/'DIMM' , 'DSMM' , 'DOMM' , 'STMM' , 'VICP' , 'EPSI' , 

+'RPM ' , 'RIMM' , 'HPMP' , 'POIS' , 'DMMM' / 

C DATA ISL,IRI/'SDEG' ,'RIMM'/,IRN/'RIMM'/ 

DATA IRN/'RIMM'/ 

PI=A.*ATAN(1. ) 

IPAR=0 

1 READ( 05, INPUTS, END=999) 

SDEG=RIMM 

IRAD=2 

IF(IWORN.NE.O)IRAD=3 

AINIT=A(NP) 

A(NP)=A(NP)-DX 

STR(8)=IRN 

C IF(IRAD.EQ.1)STR(8)=IRI 
IF(IRAD.EQ.2)STR(8)=IRN 


PLGOOOlO 

PLG00020 

PLG00030 

PLG00040 

PLG00050 

PLG00060 

PLG00070 

PLG00080 

PLG00090 

PLGOOlOO 

PLGOOllC 

PLG00120 

PLG00130 

PLG00140 

PLG00150 

PLG00160 

PLG00170 

PLG00180 

PLG00190 

PLG00200 

PLG002I0 

PLG00220 

PLG00230 
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IF(IRAD.EQ.3)STR(8)=IRN 
IFdPAR.NE.DPRINT 20,STR(NP) 

20 FORMAK / / 8X , A4 , nX , ' HMI C ' , 1 IX , ’ PMPA' , 1 IX , ' CCPM ’ / ) 

» DO 5 1=1, N 

A(NP)=A(NP)+DX 
DI=A(1)*.001 
, DS=A(2)*.001 

DO=A(3)*.001 
ST=A(4)*.001 
VI=A(5)*.001 
E=A(6)*6895. 

EN=A(7)/60. 

SL=A(8)*PI/180. 

RI=A(8)*.001 

HP=A(9)*1.E6 

PO=A(10) 

DM=A(11)*.001 

U=2.*EN*ST 

T=(DO-DI)/2. 

T2=(DM-Dl)/2. 

R=(DI+DO)/4. 

DEL=(-DI+DS)/2. 

RMX=(DI+DM)/4, 

P0=DEL*E*T/R/R+HP 
P2=DEL*E*T2 / RMX / RMX+HP 
DELO=HP*R*R/E/T 
DELT=DEL+DELO 
R2=SQRT(2.) 

H0=3.*VI*U/P0/SL 
DF=e*T’«^3 / 1 2 . / ( 1 . - PO**2 ) 

ELDF= ( R*R*DF / E / T ) ** . 2 5 
BCON=6 . *U*VI*ELDF*R*R/E/T 
BET1=BC0N/(DELT)**3 
DBAR=(RI-DEL)/DELT 
THBAR=SL*ELDF/DELT 
THEBAR=THBAR- . 200039 
THWBAR=THBAR-.5/R2 
FC0N=DF/ELDF**3*DELT 
F1NPMM=2 .4995 l*FCON*l . E-3 
F2NPMM=.42175*FCON*l.E-3 
ELF2MM=ELDF*1 . 6790*1 . E3 
SEDEG=THEBAR*DELT/ ELDF*180 . / PI 
IF(IRAD.NE.O)GO TO 222 
H 1=H0*THBAR/ THEBAR 

H0=( BCON / BETSLC BET 1 , SL , DELT / ELDF ) )** ( 1 . / 3 . ) 

PMMPA=P0*H 1 / HO* 1 . E-6 

222 IF(IRAD.EQ.l)H0=(6.*PI/16.*SQRT(2.*Rl)*VI*U/P0)**(2./3.) 

IF(IRAD.NE.2)G0 TO 90 
CALL BETCAL ( BETl , DBAR , BETA , XM , PM ) 

HO = ( BOON / BETA ) ** ( 1 . / 3 . ) 

PMMPA=PM*P0*l.E-6 

STRMPA=RI*E/R*l.E-6 

XMMM=XM*ELDF*1000. 

90 IF(IRAD.NE.3)G0 TO 60 


PLG00240 

PLG00250 

PLG00260 

PLG00270 

PLG00280 

PLG00290 

PLG00300 

PLG00310 

PLG00320 

PLG00330 

PLG00340 

PLG00350 

PLG00360 

PLG00370 

PLG00380 

PLG00390 

PLG00400 

PLG00410 

PLG00420 

PLG00430 

PLG00440 

PLG00450 

PLG00460 

PLG00470 

PLG00480 

PLG00490 

PLG00500 

PLG00510 

PLG00520 

PLC00530 

PLG00540 

PLG00550 

PLG00560 

PLG00570 

PLG00580 

PLG00590 

PLG00600 

PLG00610 

PLG00620 

PLG00630 

PLG00640 

PLG00650 

PLG00660 

PLG00670 

PLG00680 

PLG00690 

PLG00700 

PLG00710 

PLG00720 

PLG00730 

PLG00740 

PLG00750 

PLG00760 

PLG00770 
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CALL RINSBW ( ( RI +DELO ) / BELT , WP , ELCN ) 

SLW=WP*DELT/ELDF 
H0=3.*VI*U/P0/SLW 
60 CONTINUE 

IF( IRAD. EQ.4 )H0=H0*THBAR/THWBAR 
Q=U*H0* ( PI / 2 . )**3*DS / 4 . 

CCPM=Q*60.*1.E6 
HMIC=H0*1.E6 
' PMPA=P0*l.E-6 
P2MPA=P2*l.E-6 
ELDFMM=ELDF*1000. 

DELTMM=DELT*1000. 

IF(IPAR.EQ.1)G0 TO 61 
PRINT 21,A(NP),HMIC,PMPA,CCPM 
GO TO 5 

' . 61 BETA=6.*U*VI*ELDF/H0**3*R*R/E/T 

WRITE(6,1000)STR(NP),A(NP),CCPM,HMIC 
1000 F0RMAT(/1X,A4,'=' ,G15 . 7 ,5X, ' CCPM=' ,G15 . 7 ,5X, ’ HMIC=' ,G15.7) 

- WRITE (6,1 200 )PMPA,P2MPA,ELDFMM, BETA 
1200 FORMATdX, 'PMPA=' ,G15 . 7 , 5X, ' P2MPA=* ,G15 . 7 ,5X, / IX, 

+’ELDFMM=’ ,G15.7,5X, ’BETA=’ ,G15.7) 

IF(IRAD.EQ.2.0R.IRAD.EQ.3)G0 TO 888 
GO TO 887 

888 WRITE( 6 , 1 700 ) DELTMM , BETl , DBAR , STRMPA 
1700 FORMATdX, 'DELTMM=' ,G15.7 ,5X, ’BET1=' ,G15.7,5X, 

+’DBAR=' ,G15.7,5X,/1X, 'STRMPA=' ,G15.7) 

GO TO 886 

887 WRITE(6, 1400 )DELTMM, BETl, THEBAR 
1400 FORMATdX, 'DELTMM=' ,G15 .7 ,5X, 'BET1=' ,G15 . 7 , 5X , ' THEBAR= ' ,G15.7) 
WRITE(6,1600)F1NPMM,F2NPMM,ELF2MM,THBAR,THWBAR 
1600 FORMATdX, 'F1NPMM=' ,G15 .7 ,5X, 'F2NPMM=' ,G15 . 7 , 5X, / IX, 

+'ELF2MM=' ,G15.7,5X, 'THBAR=' ,G15 . 7 ,5X, / IX, ' THWBAR=' ,G15.7) 
IF(IRAD.EQ.0)WRITE(6,1500)SEDEG,PMMPA 
1500 FORMATdX, 'SEDEG=' ,G15 . 7 , 5X , ' PMMPA=' ,G15.7) 

GO TO 5 

886 IF(IRAD.EQ.2)WRITE(6,1300)XMMM,PMMPA 
1300 FORMATdX, 'XMMM=' ,G15 . 7 ,5X, ' PMMPA=' ,G15.7) 

5 CONTINUE 
21 F0RMAT(1P4E15.5) 

A(NP)=AINIT 
GO TO 1 
999 STOP 
END 

FUNCTION YLAG(XI,X,Y,IND,N1,IMAX,IEX) 

DIMENSION Xd),Y(l) 

N=N1 

IEX=0 

IF(N.LE,IMAX) GO TO 1 

N=IMAX 

IEX=N 

1 IF(IND.GT.O) GO TO 4 
DO 2 J=1,IMAX 
IF(XI-X(J))3,13,2 

2 CONTINUE 
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IEX=1 
GO TO 7 

3 IND=J 

4 CONTINUE 
IFdND.GT.DGO TO 5 
IEX=-1 

5 INL=IND-(N+l)/2 
IF(INL.GT.O) GO TO 6 
INL=1 

6 INU=INL+N-1 
IF(INU.LE.IMAX) GO TO 8 

7 INL=IMAX-N+1 
INU=IMAX 

8 S=0. 

DO 30 J=INL,INU 
IF(XI.EQ.X(J)) GO TO 13 
30 CONTINUE 
P=l. 

DO 11 J=INL,INU 
P=P*(XI-X(J)) 

D=l. 

DO 10 I=INL,INU 
IF(I.NE.J) GO TO 9 
XD=XI 
GO TO 10 

9 XD=X(J) 

10 D=D*(XD-X(I)) 

11 S=S+Y(J)/D 
YLAG=S*P 

12 RETURN 

13 YLAG=Y(J) 

GO TO 12 
END 

SUBROUTINE RINSBW(D,S,ELC) 

DIMENSION Z(33),RD(33),WP(33) 

DATA Z,RD,WP/ 

+0.0 ,0.100000,0.200000,0.300000,0.400000,0.500000,0.600000, 

+0 . 700000 , 0 . 800000 , 0 . 900000 , 1 . 000000 , 1 . 100000 , 1 . 200000 , 1 . 300000 , 
+1 . 400000 , 1 .500000 , 1 . 600000 , 1 . 700000 , 1 . 800000 , 1 . 900000 , 2 . 000000 , 
+2 . 100000 , 2 . 200000 , 2 . 300000 , 2 . 400000 , 2 . 500000 , 2.600000 , 2 . 700000 , 
+2 .800000 , 2.500000 , 3 . 000000 , 3 . 100000 , 3 . 200000 , 

+0.0 ,0.058137,0.113905,0.167663,0.219709,0.270295,0.319644, 

+0.367955,0.415412,0.462189,0.508452,0.554369,0.600107,0.645841, 
+0.691754,0.738045,0.784933,0.832661,0.881509,0.931801,0.983925, 
+1.038346,1.095644,1.156553,1.222030,1.293364,1.372356,1.461642, 
+1.565312,1.690222,1.849200,2.070620,2.439467, 

+0.0 ,0.052582,0.100776,0.145235,0.186492,0.224989,0.261100, 

+0.295148,0.327413,0.358145,0.387572,0.415899,0.443322,0.470026, 
+0.496192,0.521999,0.547635,0.573295,0.599191,0.625562,0.652681, 
+0.680876,0.710548,0.742207,0.776520,0.814399,0.857135,0.906661, 
+0 . 966037 , 1 . 040530 , 1 . 140267 , 1 .288397 ,1.557204/ 
XI=SQRT(SQRT(ABS(D-1. ))) 

IND=0 

S=YLAG(XI,RD,WP,IND,3,33,IEX)**3 
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ELC=YLAG(XI,RD,Z,IND,3,33,IEX) 

RETURN 

END 

SUBROUTINE BETCAL(BET1 , DEAR, BETA, XM, PM) 
DIMENSION HL1(13),CL(13) 

DIMENSION HL(7),WK1(30),HLP(1),BLP(1),ANS(1,7) 
DIMENSION BQ(8),AQ(8,7) 

DIMENSION BL(1 1 ) , PMAX( 11,7) ,XPM( 11,7) 

DATA HL1,CL/ 


DATA BQ,AQ/ 

+3.68888,4.60517,5.29832,6.21461,6.90776,7.60090, 
+1.3849,1.1936,1.1249,1.0732,1.0491,1.0314,1.0152 
+1.2394,1.1603,1.0927,1.0590,1.0355,1.0149,1.0055 
+1.1937,1.1080,1.0660,1.0373,1.0134,1.0032,1.6240 
+ 1.1257, 1.0743 , 1.0395 ,1.0120, 1.0009 ,1 . 7075 , 1.4036 
+1.0801,1.0418,1.0117,0.9999,1.7801,1.4427,1.2836 
+1.0448,1.0123,0.9996,1.8471,1.4800,1.3078,1.1625 
+1.0137,1.0000/ 

DATA BL/ 

+3.91202, 4.60517, 5.29832, 6.21461, 6.90776, 7.60090, 8.51719, 
+ 9.21034, 9.90349,10.81978,11.51293/ 

DATA PM AX/ 


HLP(l)=ALOG(DBAR) 

IND=0 

CON=EXP( YLAG ( HLP ( 1 ) , HLl , CL , I ND , 3 , 1 3 , I EX ) ) 
B ET AC= ( CON / BET 1 ) ** ( 1 . / . 8 1 9 5 ) 

BETA=BETAC 
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+ 0.33611, 
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0.91745, 
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1.99386, 
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+ 3.06409, 

3.60775, 

4.30551, 

4.85468, 

5.70445, 

6.89099/ 
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,1.0068,1.4179, 

,1.5028,1.2931, 

,1.3592,1.2331, 

,1.2599,1.1381, 

,1.1500,1.0865, 

,1.0938,1.0488, 
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IF(BETA.GT.1.E4)G0 TO. 5 
BLP(l)=ALOG(BETAC) 

CALL IBCIEU(AQ,8,BQ,8,HL,7,BLP,L,HLP,1,ANS,1,WK1,IER) 
BETA=ANS(1,1)*BETAC 
5 BLP(l)=ALOG(BETA) 

CALL IBCIEU(PMAX,11,BL,11,HL,7,BLP,1,HLP,1,ANS,1,WK1,IER) 
PM=ANS(1,1) 

CALL IBCIEU(XPM,11,BL,11,HL,7,BLP,1,HLP,1,ANS,1,WK1,IER) 
XM=ANS(1,1) 

RETURN 

END 

FUNCTION BETSL(BETl,SL,DOL) 

DIMENSION BL(13),BFN(13) 

DATA BL/ 

+13.12236,12.20607,11.51293,10.81978, 9.90349, 9.21034, 8.51719, 
+ 7.60090, 6.90776, 6.21461, 5.29832, 4.60517, 3.91202/ 

DATA BFN/ 

+-2.42400,-2.27362,-2.15682,-2.03710,-1.87397,-1.74656,-1.61545, 

+-1.43609,-1.29564,-1.15098,-0.95330,-0.79907,-0.64119/ 

PAR=ALOG ( . 5*BET1**( 2 . / 3 . )*DOL/ ( SL-DOL/SQRT( 2 . ) ) ) 

IND=0 

BETSL=EXP(YLAG(PAR,BFN,BL,IND,3,13,IEX)) 

RETURN 

END 
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